Astrophysics and Space Science 
DOI 10.1007 /$•••••-•••-••••-• 



The effect of two-temperature post-shock accretion flow 
on the linear polarization pulse in magnetic cataclysmic 
variables 



Gordon E. Sarty • Curtis J. Saxton • Kinwah Wu 



OO 

O 

o 

<: 
in : 



^ . 

Or 
i < 

O ■ 

to: 



> 

m 

in 
q 

od 
o 

OO 

o 



© Springer- Verlag •••• 

Abstract The temperatures of electrons and ions in 
the post-shock accretion region of a magnetic cata- 
clysmic variable (mCV) will be equal at sufficiently 
high mass flow rates or for sufficiently weak magnetic 
fields. At lower mass flow rates or in stronger magnetic 
fields, efficient cyclotron cooling will cool the electrons 
faster than the electrons can cool the ions and a two- 
temperature flow will result. Here we investigate the 
differences in polarized radiation expected from mCV 
post-shock accretion columns modeled with one- and 
two-temperature hydrodynamics. In an mCV model 
with one accretion region, a magnetic field <L 30 MG 
and a specific mass flow rate of ~0.5 g cm -2 s _1 , along 
with a relatively generic geometric orientation of the 
system, wc find that in the ultraviolet either a single 
linear polarization pulse per binary orbit or two pulses 
per binary orbit can be expected, depending on the ac- 
cretion column hydrodynamic structure (one- or two- 
temperature) modeled. Under conditions where the 
physical flow is two-temperature, one pulse per orbit 
is predicted from a single accretion region where a one- 
temperature model predicts two pulses. The intensity 
light curves show similar pulse behavior but there is 
very little difference between the circular polarization 
predictions of one- and two-temperature models. Such 
discrepancies indicate that it is important to model 
some aspect of two-temperature flow in indirect imag- 
ing procedures, like Stokes imaging, especially at the 
edges of extended accretion regions, were the specific 
mass flow is low, and especially for ultraviolet data. 
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1 Introduction 



polarization • white 



Magnetic cataclysmic variables (mCVs), consisting of 
the classes known as polars and intermediate polars, 
are composed of a Roche-lobe-filling M type main se- 
quence s tar in orbit about a magnetic white dwarf (see 
Warner ( 19951 ) for a review). Mass is lost through the 
inner Lagrangian point, L\, and flows toward the mag- 
netosphere of the white dwarf either predominately in a 
stream (polars) or after forming a truncated accretion 
disk circulating around the white dwarf (intermediate 
polars). In either case, the ionized gas follows mag- 
netic field lines to the surface of the white dwarf after 
the gas reaches the magnetosphcrc where the magnetic 
pressure exceeds the gas ram pressure. Upon reach- 
ing the white dwarf surface the gas will be essentially 
at "free fall", with highly supersonic velocities. The 
abrupt stop of the radial inflow near the surface of the 
white dwarf leads to the formation of a shock, which 



heats the inflowing material ( 


Fabian. Prinele and Reed 


1976t Kine and Lasota 


19791: Lamb and Masters! Il979t 


Wull2000l). The hot subsonic post-shock flow settles 



gradually onto the white dwarf, and cools via emitting 
brcmsstrahhmg X-rays and optical/infra-red cyclotron 
radiation. 

The hydrodynamic structure of the post-shock set- 
tling flow is determined by radiative and particle en- 
ergy processes, which are essentially characterized by 
the bremsstrahlung cooling time tfa, the cyclotron 
cooling time t C y, the electron-ion energy-exchange 
time t P \ , the electron-electron collisi onal time tee. 



and t he ion-ion collisional time t\\ ( King and Lasota 



1979t 



Saxton et al 



Lamb and Masters! Il979t Imamua et al. 1996 : 
For weakly magnetic systems 



2 



(with B ~ 10 6 G or weaker) 
erg cm~ 2 



nositics 10 31 - 10 33 



tcy > ^br ^ ^ei ^ * 



with accretion lumi- 
s _1 , typical of mCVs, 
As the strength of the magnetic 
field increases to B ^ 10 7 G, cyclotron cooling may 
dominate bremsstrahhmg cooling, t C y < t^ r . For suffi- 
ciently strong magnetic fields and low specific accretion 
rates, i C y is so short (t C y < t e {) that collisions between 
electrons and ions cannot maintain an equal tempera- 
ture between the two types of particles. The accretion 
flow is therefore in a two-fluid regime which requires 
a two-temperature (2T) hydrodynamic description. A 
strong magnetic field can also result in a situation where 
cyclotron radiative loss is so rapid (t C y < tee) that 



cy 

electron-electron collisions are not efficient enough to 
maintain a Maxwellian distribution. In the extreme sit- 
uation where t C y < ^ib the accretion flow is no longer 
hydrodynamic . 

Previous calculations of cyclotron radiation from the 
post-shock settling flow in mCVs have either assumed a 
uniform density and temperature ( Chanmueam and Dulk^ 
( 198ll); [Meggitt and Wickramasinghel dl982l); 

Barrett and Chanmugam ( f984 ); 

Wickramasinghe and Meggittl (If 9851) ; 

Wu and Chanmugaml (If 988l.lf 9891) ) or a one-temperat ure! 

(IT) structu r e rtWu and Wickramasinghe! (|l990Lll992h : 



radiative transfer calculations, and in |J4] we examine 
the differences between the spectro-polarization prop- 
erties of the emission from IT and 2T accretion flows 
and discuss their implications. Concluding remarks are 
made in 



2 Method 

Calculations were done for 36 cases: 3 values of mag- 
netic field strength (B-j = f, 3 and 5) x 3 values of 
white dwarf mass (Af W D = 0.5, 0.7 and f.O M ) x 2 
values of mass flow rate (M = fO 16 and f0 14 g s _1 ) x 2 
values of column hydrodynamic structure (IT and 2T). 
Here Br is the magnetic field strength at the accretion 
spot in units of f 7 G. In all cases the area of the accre- 
tion spot, A, was taken to be constant at 2 x 10 14 cm 2 . 
This value for the area is based on the consideration 
that 



Potter et all (|2002h ). However, detailed ID calcula- 
tions of the hydrodynamic structure of a post-shock 
accretion column that self-consistently include cy- 
clotron and bremsstrahlung cooling clearly show that 
a 2T structure is to be expe cted in many p hysica l 
situations relevant to mCVs (llmamua et al.l ([1996) ; 



20011) ; 



Woelk and Beuermannl ( 19961); ISaxt on and Wu 
Wu et all (|2003l ): ISaxton et all (|2005l )V 

Here we have computed and compared the cyclotron 
radiation from a cylindrical post-shock accretion col- 
umn, with a uniform cylindrical radial structure, as- 
suming both a IT hydrodynamic structure and a 2T 
hydrodynamic structure. The resulting cyclotron spec- 
tra for a grid of three white dwarf masses (0.5, 0.7 and 
1.0 Mq), three magnetic field strengths (10, 30 and 50 
MG) and two mass flow rates (10 16 and 10 14 g s _1 ) 
were computed for various viewing inclination angles. 
For each case, using the computed viewing- angle de- 
pendent cyclotron spectra, Johnson bandpass (jjohnson 
196.4 lBessel3ll990h filtered light curves over an orbital 
period were computed for a mCV with an orbital incli- 
nation of 45° and with the given accretion column at a 
co-latitude of 30°. 

This work is organized as follows. In ^ we out- 
line the hydrodynamic formulation used to determine 
the density and temperature structure of the post-shock 
flows and the radiative transfer through the ionized ac- 
creting gas. In $3]we present the results of the polarized 



A = 4^< D x (1) 

where i? W D is the white dwarf radius and x ~ 10 was 
chosen as being a typical value for the fraction of white 
dwarf surface occupied by the accreti on spot. The whit e 
dwarf radius, in turn, is taken to be (|Nauenberg||l972l) : 



i?wo _ 0.0225 s/i - (M wv /M 3 )*/3 



R, 



Mwd (M WD /M 3 ) 1 /3 



(2) 



where M3 = 5.816M //i 2 VD is the Chandrasekhar mass 
limit and /j WD is the mean molecular weight of the white 
dwarf material taken to be equal to 2.00. The chosen 
spot area, A, gave specific accretion rates, m = M/A, 
of 50.0 g cm- 2 s- 1 for M = 10 16 and 0.5 g cm^s" 1 
for M = 10 14 . With this constraint on the specific 
mass accretion rate, IT and 2T hydrodynamic accretion 
column structures were computed as described in i j2.ll 
The computed electron temperature and density in the 
column, as a function of height, was then determined 
for the center of each cube, of individual side length 
Ah, in a computational grid-of-cubes. Each individual 
cube was modeled as having a uniform temperature and 
density at the center values. Ray-tracing was then done 
through the computational grid-of-cubes as described in 
32.21 The radiation transfer calculation method for the 
polarized cyclotron radiation from the accretion column 
emerging from the top of the grid-of-cubes for various 
column orientations (viewing angles) is also described 
in 32.21 The method for computing a broadband light 
curve over the course of a binary orbit is given in 



mCV two-temperature flows 
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2.1 Hydrodynamics 

The hydrodynamic structure of the post-shock accre- 
tio n column was c ompu ted following the methods given 



by ISaxton et alj ((2005) . The hydrodynamic structure 
for the two-temperature flow is specified by the conser- 
vation of mass, momentum, energy of the electrons and 
energy of the ions, respectively, as the height above the 
white dwarf, z, changes according to the relationships 
explicitly given by 



dp du 
u-f + p—=0 

dz dz 



du 
pu— = 

dz 



dP e dPj 

dz dz 

dP e uP e dp 
u— 7 — = (7 

dz p dz 



i)(r ei -A) 



dp { 

dz 



7 



uP[ dp 
p dz 



-(7-i)r e 



(3) 



(4) 



(5) 



(6) 



where u is the flow velocity, p is the density, P e is the 
electron partial pressure, Pj is the ion partial pressure, 
F e j is the electron-ion energy exchange rate, A is the 
electron cooling rate and 7 is the adiabatic index. An 
ideal gas law for mono-atomic species was assumed for 
both electrons and ions so that 7 = 5/3. The total 
pressure is P — P + Pj and the number densities, rij, 
are given by Pj = rijkTj where k is the Boltzmann 
constant and j — c or i. The cooling rate is the sum of 
the cooling due to bremsstrahlung and cyclotron emis- 
sions, A = [A br + A cy ] oc [(1/tbr) + (l/*cy)]- In rel- 
ative terms, the cyclotron cooling is most efficient in 
the hot, less dense plasma just below the shock and the 
bremsstrahlung is more efficient at the base of the accre- 
tion column where the the gas is dense. The electron- 
ion energy exchange rate T e j oc 1 /t e j is a function of the 
electron and ion number densities and temperatures. 
One-temperature flow is similarly given by 



dp du 
u-f + p—=0 

dz dz 



dP 

dz 



du 
pu— = 

dz 



dP uP dp , . 
u 7 — = _( 7 _ i)A. 

dz p dz 



(7) 



(8) 



(9) 



The flow is assumed to have a constant density, 
pressure and temperature across the column area 



A, Note that variation of structure across the col- 
umn can substantially affect the rad iation produced 
jWu and Wickramasinghlll990l . fl992l) but here we as- 
sume the simpler structure so that the effects of the 
IT and 2T structures, alone, in the production of ra- 
di ation may be discer ned. In the formulation give n 
by Saxton et al. ( 2005[ ) (see also Imamua et al. ( 19961 )) 
the pressure ratio of the electrons to ions at the shock, 
o"s = Pe,s/Pi s was taken to be a free parameter. Here 
we take a s — Zq = 1.099, which corresponds to an 
assumption of solar metallicity. 

The results of the hydrodynamic calculations for 
the column profile were transferred to a computational 
grid-of-cubes of size Nt X N X N as the column was ro- 
tated in that grid. More exactly, the temperatures and 
densities of a radially uniform cylinder were computed 
at the center of each cube in the grid and those values 
were taken to be uniform over the individual cube. The 
distance between cube centers on the grid, Ah, was uni- 
form in the x c , y c and z c grid directions and was deter- 
mined by dividing the computed column (shock) height, 
h s , by a fixed number of strata so that the spatial res- 
olution in the z direction of the column was fixed. The 
number of strata was set at 20 for all cases except in six 
(unimportant, it turned out) cases where it resulted in 
a grid size so large that subsequent computations were 
impractically long. Those six cases were the high mass 
flow rate (M = 10 16 g s _1 ), M WD = 0.5M Q cases where 
5 strata were used and where it turned out that there 
were no differences between the IT and 2T cases. With 
the number of strata fixed, the computational grid-of- 
cubes size was determined as follows. The extent of the 
grid in the y c and z c directions was taken as \c\ (mean- 
ing round c up to the nearest full multiple of Ah) where 
c = y/hg + d 2 and d = I^JApK where A is the column 
spot area of 2 x 10 14 cm 2 . The thickness of the grid 
in the x c direction was \d~\ . With the grid sized this 
way, the accretion column could be rotated completely 
within the grid and the ray-tracing through the column 
was always done vertically through the grid, in the z c 
direction, with the computed radiation emerging from 
the top x c y c surface of the computational grid-of-cubes 
as shown in Fig. [T] The resulting grid size, in number 
of grid cubes, was then N t x N x N . With the shorter 
columns this size could be quite large, on the order of 
300 x 300 x 300. 

2.2 Polarized radiative transfer 

After the values of electron temperature and number 
density, T e (z) and n e (z), were determined for the z 
(height above the white dwarf surface) value at the cen- 
ter of each cube in the computational grid, the optical 



1 




Fig. 1 The orientation of the post-shock accretion column in the ray-tracing computational grid of cubes for a tall and 
short column for various viewing angles (indicated above the images). The greyscale within the columns is proportional to 
the computed electron temperature with white being hot. (The greyscale is relative to the hottest temperature for each 
column; the temperatures of the two different columns cannot be compared directly to each other simply by looking at 
the greyscale.) The side view, in section through the middle, of the columns are shown and the columns are circular in 
cross-section. Both columns are from the Mwd = I.OMq cases. The tall column, shown in panels (a)-(c), is from the 2T 
calculation with B7 = 1.0 and M = 10 14 g s _1 (m = 0.5 g cm~ 2 s -1 ). The short column, shown in panels (d)-(f), is from 
the 2T calculation with B 7 = 1.0 and M = 10 16 g s _1 (m = 50 g cm^s -1 ). The height of the tall column is 30.2 x 10 6 
cm and the height of the short column is 3.2 x 10 6 cm; both columns have a spot area A — 0.20 x 10 15 cm 2 . Both columns 
are divided into 20 strata along the z direction with a resulting grid size N t x N x TV = 13 x 29 x 29 for the tall column 
and Nt x TV x TV — 101 x 107 x 107 for the short column. Both grids have the same number of temperature bins in the 
z direction but the aspect ratio of the short column required a much larger grid, and corresponding computation time, 
for the ray tracing radiative transfer calculations. The surrounding black boxes in the images represent the ray-tracing 
computational grid of cubes (the individual cubes are the size of pixels in these images) and ray tracing was always done 
in the z c direction; the computational box was always viewed from above when summing the emission from the individual 
cubes. 
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depth parameter S for each cube was computed accord- 
ing to 



S(z) = 2.01x10 s 



Ah 
10 5 cn 



n e (z) 



3 x 10 7 G 



B 



(10) 



where B is the magnetic field streng th 
(jWickramasinghe and Meggittll985h . finally rays were 
traced from the bottom to the top of each (x c , y c ) col- 
umn of the computational grid using a finite difference 
approximation (i.e. dz c — Ah) of 
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(11) 



( Pacholczvk and Swihartj 19751 ). where I v , Q U: U v and 
V v are the (frequency dependent) Stokes intensities (see 
e.g. Rvbicki and Lightman ( 19791 )). k v , q v , and v v are 
the respective corresponding absorptivities, f v and h v 
are the Faraday mixing coefficients and e/.„, cq,^, eu,v 
and ev> are the corresponding emissivities. The ab- 
sorptivities and emissivities are related by Kirchhoff's 
law as the accretion column's electrons are in local 
thermodynamic equilibrium and the absorptivities and 
e missivities are (essentially) functions of 5, T ei v and 
8 (jMeggitt and Wickramasinghelll982l ). where 8 is the 
angle between the magnetic field direction (z axis) and 
the viewing direction [z c axis). The Stokes intensities 
emerging from the top of each column of cubes were 
then summed over the columns to give the four total 
intensities emerging from the top of the grid-of-cubes. 

From the Stokes intensities, the linear polarization 
L v = y Q u + U% jl v and circular polarization C„ = 
V v /I v were computed for 18 angles between 0° and 
90° and for 25 frequencies, expressed in cyclotron har- 
monic number H = lu/lo c , satisfying 2 < H < 12. Here 
lo c = eB/m e c is the cyclotron frequency where e and 
m e are the electron charge and mass, respectively, B is 
the magnetic field strength and c is the speed of light. 

The ray tracing calculations were done on a Beowulf 
cluster of 20 computers with each node in the cluster 
working separately on one of the 30 cases with 20 cylin- 
der strata. Each node consisted of a 1 GHz Pentium 
III processor with 512 MB of RAM. The smaller grids 
required under a day to complete while the larger grids 
took several weeks of CPU time. The 6 large problems 
where 5 cylinder stata were used were parallelized to 
compute one viewing angle per node to reduce com- 
puting time and even then it took about 4 weeks to 
complete the computations. The code was optimized 



to minimize the use of computer memory, so that the 
entire computational cube's data did not have to be 
stored, but no directed effort was given to optimize the 
speed of the computations; an improvement in the in- 
terpolation of the column temperatures and densities 
to the cube centers would likely increase the computa- 
tional speed significantly. 

2.3 Orbital polarization light curves 

To better assess the difference in observational predic- 
tions between a IT and a 2T accretion column, the 
intensity, linear and circular polarization of the emis- 
sion from the column as seen from Earth over a binary 
orbital period was computed. For this the angle, 8, be- 
tween the viewing direction and the magnetic field at 
each angular position of the white dwarf relative to the 
mass donor star, was needed. Taking the magnetic field 
to be normal to the white dwarf surface at the accretion 
column, the viewing inclination angle, 8, is given by 



sin i cos a sin < 



(12) 



where i is the inclination of the binary orbit and <f> is 
the co-latitude of the accretion spot on the white dwarf 
surface. The orbital position is a with a = represent- 
ing inferior conjunction of the mass donor star (mass 
donor star between the distant observer and the white 
dwarf) when the accretion spot is located directly below 
the white dwarf rotation axis as seen by the observer. 
We assume that the rotation of the white dwarf is syn- 
chronized with the orbital period as this is on e of the 
defin ing features for CVs classified as polars (jWarner 
19951 ). 

The angle of linear polarization, \t, was also com- 
puted over a binary orbital period. The angle \t = 
Xv + 8 where Xv is the linear polarization angle relative 
to the computational grid-of-cubes and 6 is the angle 
that the grid-of-cubes needs to be rotated by, about the 
z c axis, to line up with the projection of the accretion 
column on the sky. Taking the orbital axis as defining 
zero position angle we have 



tan2Xt/ = U V /Q V 
and 

g 

tan(<5 -it/2) = - 



cos i cos a sin ( 



sm a smt 



(13) 



(14) 



A nominal configuration of <f> = 30° and i = 45° was 
chosen to illustrate a relatively wide range of viewing 
angles without the accretion spot disappearing behind 
the limb of white dwarf. 

The cyclotron harmonic range between 2 and 12 cor- 
responds to a wavelength range of 5350 to 892 A for 



6 



B r = 1, 1785 to 298 A for B 7 = 3 and, 1070 to 178 A for 
B 7 = 5. With the exception of the higher frequencies 
in the B 7 = 5 case, these wavelength ranges are cov- 
ered by the LKJIRVBU filte rs of the Johnson system 
I Johnsonl Il965l iBesselll ll99oT) . The response of those 
filters, in cases where the IT and 2T solutions were dif- 
ferent, was multiplied against the computed spectra to 
give Stokes vector bandpass intensities for each point 
in the orbital light curve. 



3 Results 

In all cases, substantial linear polarization was evident 
generally only at high viewing angles, where the accre- 
tion column is viewed side-on. so it was at high viewing 
angles where the largest differences between the linearly 
polarized emission from IT and 2T columns would be 
seen. Cyclotron harmonic spike-features were visible in 
the linear polarization spectra at low harmonic num- 
bers, where the emission is marginally optically thick 
at the spectral harmonic peaks and is optically thin in 
between the spectral harmonic peaks. The linear po- 
larization spectra became smooth at higher harmonic 
numbers and the cyclotron harmonic spike-features dis- 
appear, since the emission becomes optically thin at 
shorter wavelengths. At the high accretion rate of 
M = 10 16 g s _1 , there was no obvious difference in 
the linear polarization spectra of the IT and 2T models. 
This high accretion rate is sufficient for the electron-ion 
collisions to maintain approximately equal ion and elec- 
tron temperatures, in spite of radiative loss. Circular 
polarization was most prominent when viewed along the 
accretion column (magnetic field) at low viewing angles, 
so it was at low viewing angles where the largest dif- 
ferences between the circularly polarized emission from 
IT and 2T columns would be seen. The relative differ- 
ence between the predicted circularly polarized radia- 
tion by the IT and 2T hydrodynamic formulations was 
negligible, as the accretion flow is essentially IT given 
the efficient energy exchange between the electrons and 
the ions. At the lower mass flow rate of M = 10 14 g 
s _1 , the situation is very different. Electron-ion col- 
lisions become less efficient at lower densities. For a 
sufficiently strong magnetic field, electrons will lose en- 
ergy rapidly via cyclotron radiation but their collisions 
with ions are unable to maintain an equal temperature 
between the two particle species. For B 7 = 3 and 5, 
the flow is essentially 2T, and we find very different 
predicted polarization for the IT and 2T formulations. 

To investigate the observational difference in the pre- 
dictions of a IT and 2T accretion column for orbital 



light curves, the Johnson filter with the highest fre- 
quency pass-band that contained the computed emis- 
sion, where the IT and 2T emission differences were 
the greatest, was chosen for the B 7 = 1 and 3 cases. 
For Br = 1 the filter choice corresponds to the infra- 
red J band, for B 7 = 3 the choice corresponds to the U 
band. For the B 7 = 5 case, the light curve at the fre- 
quency corresponding to the cyclotron harmonic 9.92 
(the 20th of the 25 frequencies computed), in the ul- 
traviolet beyond the U filter, was computed. These 
choices emphasize the differences between the IT and 
2T model predictions at the higher cyclotron harmon- 
ics. Since polars are not generally observed in the ultra- 
violet, V band orbital light curves for the B 7 — 3 and 5 
cases, where the flow is essentially a 2T flow, were also 
computed. 

With the given choice of filters, the differences in lin- 
ear polarization, linear polarization angle, intensity and 
circular polarization light curves are shown in Figs. [2] 
to [5] The IT and 2T predicted polarization angle 
and circular polarization light curves are generally the 
same, but there are some differences for the high white 
dwarf mass, strong magnetic field cases in the ultra- 
violet. There arc significant differences between the 
model predictions for the linearly polarized and inten- 
sity light curves in the ultraviolet. In two cases the 
IT model predicts two pulses per orbit where the 2T 
model predicts one pulse, for both linear polarization 
and intensity, with an increase in linear polarization 
and intensity at the time of superior conjunction when 
we see the accretion column at the top of the white 
dwarf as seen from the observer. The two cases are: 
B 7 = 3,M WD = 0.7M Q and B 7 = 5, M WD = 1.0M Q . 
However, in the visual V band there are no appreciable 
differences in the orbital light curves predicted by the 
IT and 2T models. 



4 Discussion 

The emission of optical/infra-red polarized light from 
magnetic CVs is generally thought to arise from the 
post-shock region of the magnetically funneled accre- 
tion flow onto the surface of the white dwarf. How- 
ever the structure of that accretion column may be 
considerably more complex than the simple cylin- 
der considered here. For example, it is probable 
that the threading of the gas stream from the L\ 
point onto the magnetic field lines at the magneto- 
sphere boundary happens over a range of orbital lon- 
gitudes. This spread of the threading region leads to 
a spr ead in the footprint of the accretion column into 
arcs ( Ferrario and Wickramasinghel fl99oh or possibly 
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J Filter 
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(m) p"" 6 (n) p"" 6 (o) i"" 15 ' 

Fig. 2 Linear polarization orbital light curves comparing predictions for the one-temperature model ("IT", red line) with 
the two-temperature ("2T", green line) model. Only the low mass flow rate cases, M = 10 14 g s" 1 (to = 0.5 g cm 2 s 1 ), are 
shown as there are no differences in the two model predictions at the higher mass flow rate. Light curves through bandpass 
filters at the higher harmonics are shown in the first nine panels (a) to (i). The cases in which the 2T model produced 
different radiation output predictions from the IT model are shown again in the last six panels (j) to (o) for the visual V 
filter. 
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Fig. 3 Linear polarization angle orbital light curves comparing predictions for the one-temperature model ("IT", red 
with the two-temperature ("2T", green line) model. The figure is laid out in the same way as Fig. [2] 
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Fig. 4 Intensity orbital light curves comparing predictions for the one-temperature model ("IT", red line) with the 
two-temperature ("2T", green line) model. The figure is laid out in the same way as Fig. 
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Fig. 5 Circular polarization orbital light curves comparing predictions for the one-temperature model ("IT", red 
with the two-temperature ("2T", green line) model. The figure is laid out in the same way as Fig. [2] 
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lex patterns 


(Potter. Hakala and Cropper 


Potter et al. 


2004) on the surface of the 



white dwarf. The complexity of the accretion col- 
umn has undoubtedly contributed to the difficultly 
in modeling a g ood fit to obser v ed polarimetr i c ligh t 



curve data (e.g. ICropper] ill 9891): IPiirola et all (Il990h ; 
Wickramasinghe et all |l99ll ): IPotter et all (jl997h 
Buckley et al I (|2000l )). especially to linear polarization 
data. Here we have shown that the detail of model- 
ing the post-shock flow as a IT or 2T flow can have a 
substantial effect on the predicted polarization, partic- 
ularly at lower mass flow r ates and at m agnetic fields 
strengths typical for polars (|Warnerlll995l ) . This differ- 
ence in predicted radiation would matter, in particular, 
for the lower specific mass flow rates required in models 
of extended accretion regions, especially at the edges of 
the accretion regions. 

For the lower mass flow rate of 0.5 g cm _2 s~ 1 , we 
found that the flow was essentially IT for B7 = 1 and 
2T for the stronger fields, B7 = 3 and 5. This may 
be understood on the basis of the relative magnitudes 
of the cyclotron cooling times and the electron-ion en- 
ergy exchange times £ C y and i e j ( Lamb and Masters] 
19791) . Two-temperature flows are expected to occur 
for a larger magnetic field given a fixed mass flow rate 
or for smaller mass accretion rates given a fixed mag- 
netic field because the cyclotron cooling efficiency in- 
creases relative to the electron-ion energy exchange ef- 
ficiency in those circumstances. In line with those ex- 
pectations, larger differences in the predicted orbitally- 
phascd light curves, between the IT and 2T cases, are 
seen at stronger magnetic fields where the overall cy- 
clotron cooling is more efficient. 

For differences in predicted linear polarization pulses, 
particularly note the cases of Fig. [5] (e) and (i), which 
show one predicted linear polarization pulse per orbit 
for the "correct" 2T flow versus two predicted linear 
polarization pulses per orbit when an inappropriate IT 
flow is assumed. In that case one would incorrectly 
interpret an observation of two pulses per orbit as com- 
ing from one accretion column, caused by viewing the 
accretion column sticking out from the white dwarf 
limb, if a IT flow was used to model the data. In 
actuality, for a generic white dwarf of mass 0.7 Mq 
and magnetic field of 30 MG, two linear polarization 
pulses could be produced by two accretion columns, 
one at each magnetic pole of the white dwarf. This 
consideration is important in indirect imaging cases 
where one is trying to model the white dwarf surface 
emission and determine the magnetic field configuration 
through a forward-modeling approach to the inversion 
of polariz ed light curve data in a technique l ike 'Stokes 
imaging' (jPotter. Hakala and Cropperlll9 98). To date 



Stokes imaging has assumed a homogeneous tempera- 
ture and density over the height of the column but we 
have shown that the assumption of a IT or 2T flow in 
the column can significantly alter the amount linearly 
polarized light that a forward model would give, es- 
pecially at shorter wavelengths. More accurate Stokes 
imaging models may need to incorporate, somehow, the 
fact that the flow will be 2T, especially at the rim of an 
accretio n spot, where the specific accretion rates may 
be low (jAchilleos. Wickramasinghe and Wulll992l ). 

On the positive side, for Stokes imaging, is that 
the differences in emission predicted from IT and 2T 
columns is insignificant in all cases for lower cyclotron 
harmonics in the infra-red and optical. So for infra-red 
and optical frequency light curves, as have been used 
for past applications of Stokes imaging, it may not be 
necessary to model a two-temperature flow. However 
we should note that our grid of models does not cover 
field strengths between B7 = 1 and = 3 where there 
may be small 2T effects at optical wavelengths. Also 
the variation of the accretion spot area, A, to larger val- 
ues would increase 2T effects by reducing the specific 
accretion rate. In any case, the addition of ultraviolet 
light curve data to optical light curve data for model- 
ing the emission from polars may give a truer image 
of the white dwarf surface. Such polarized UV data 
may become available from future space observatories 
like the proposed World Space Observatory for Ultra- 



violet (WSO-UV) (jPagano et all 120071 : lUslenghi et al 
20081) . We note, too, that there are essentially no dif- 



ferences in the prediction of circular polarization light 
curves between IT and 2T models. Only linear po- 
larization differences are apparent. Linear polariza- 
tion light curves were more sensitive to assumptions 
made about the accretion flow than circular polariza- 
tion light curves. This distinction between predicted 
linear and circularly polarized emission is consistent 
with past modeling of observational data, which either 
explicitly or implicitly assume one-temperature hydro- 
dynamics, where the fits to linear polarization data are 
generally p oorer than fits to the circular p ol arization 
data (e .g. Barrett and Chanmugam J 19841): ICropper 
( 1989h : IPotter et all l|l997t ): iBucklev et all ()2000h T 



5 Conclusion 

Assuming a one-temperature hydrodynamic post-shock 
accretion column as the source for polarized radiation 
in models of magnetic cataclysmic variables can lead 
to erroneous predictions for the radiation when the cy- 
clotron cooling efficiency is greater than the electron- 
ion energy exchange efficiency. This effect shows up 
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at the lower mass flow rate modeled here (to = 0.5 g 
cm _2 s _1 ) at higher cyclotron harmonics for the fairly 
generic white dwarf masses of 0.7 and 1.0 M Q and mag- 
netic fields of 30 MG or greater. So the interpretation 
of light curve data obtained at higher frequencies, in 
the ultraviolet, needs to take into account the effect of 
two-temperature flow on the production of polarized ra- 
diation. In particular, the number of linear polarization 
pulses observed in the ultraviolet can be misinterpreted 
if a one-temperature accretion column flow is assumed. 
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